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ABSTRACT 



r~| . Stellar spectroscopy, UBV photometry, Ha imaging, and analysis of data from the ATCA 

Qh' H I survey of the LMC are combined in a study of the LMC OB association LH 72 and its 

surroundings. LH 72 lies on the rim of a previously identified H I shell, SGS-14, and in the 
$H ' interior of LMC-4, one of the LMC's largest known supergiant shells. Our analysis of the H I 

C/3 . data finds that SGS-14 is expanding with velocity Vexp ~ 15 km s~^, giving it an expansion age 

' of ~15 Myr. Through the stellar spectroscopy and photometry, we find similar ages for the oldest 

stars of LH 72, ~15— 30 Myr. We confirm that LH 72 contains an age spread of ~15— 30 Myr, 
similar to the range in ages of stars derived for the entire surrounding supergiant shell. Combining 
analysis of the O and B stars with Ha imaging of the H II region DEM 228, we find that DEM 
' 228 accounts for only 60% of the available ionizing Lyman continuum photons. Comparing the 

distribution of ionized gas with that of the H I, we find that DEM 228 and LH 72 are offset by 
~ 1 — 2' from the peak 21-cm emission, towards the interior of SGS-14. Taken together, these 
results imply that SGS-14 has cleared its interior of gas and triggered the formation of LH 72. 
On the basis of our results, we suggest that LMC-4 was not formed as unit but by overlapping 
shells such as SGS-14, and that LH 72 will evolve to produce a stellar arc similar to others seen 
within LMC-4. 

Subject headings: Magellanic Clouds — ISM: bubbles — H II regions — ISM: HI — stars: early-type 

1. Introduction most striking feature is a kpc-sized ring of H II 

regions, dubbed LMC-4 by Meaburn (1980). The 
The Constellation HI region in the LMC has ^^^^ of LMC-4's interior is almost empty of other 
drawn considerable interest over the last 50 years, ionized gas, except for a single bright H II region, 
for a number of reasons. One is that the region ^-^^ "eighth-note nebula" DEM 228 (Davies, El- 
contains several structures noteworthy enough to jio^t, & Meaburn 1976). 21-cm maps, such as the 
be given names. At the wavelength of Ha, the high-resolution map of Kim et al. (1998), show 
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that LMC-4 is enclosed by a near-circular ring of 
H I gas, named SGS-U by Kim et al. (1999). The 
ring shape of H I is only spoiled by a tongue of gas 
that joins DEM 228. Lucke & Hodge (1970), in 
their mapping and analysis of OB associations in 
the LMC, identified many associations contained 
within the boundaries LMC-4. Some of these are 
organized into circular arcs of stars, two of the 
most prominent of which were named "Quadrant" 
and "Sextant" by Efrcmov & Elmegreen (1998). 
The colorful name Constellation III, as originally 
used by McKibben Nail & Shapley (1953), des- 
ignates the Sextant arc (Efremov & Elmegreen 
1998). 

More important than the names are the variety 
of physical mechanisms suggested to have created 
the kpc-sized hole in the interstellar medium and 
the substructures within (see review in Efremov 
& Elmegreen (1998)). The size of the hole sug- 
gests a large input of energy required to clear it, 
perhaps produced in a "super-supernova" (West- 
erlund & Mathewson 1966), a gamma ray burst 
(Efremov, Elmegreen, & Hodge 1998), or a high- 
velocity cloud impact with the LMC disk (Braun 
1996). Dopita, Mathewson, & Ford (1985) of- 
fered a less dramatic solution, in which star for- 
mation has propagated outwards from a central 
location, clearing a cavity through the combined 
action of supernovae and stellar winds. The pre- 
dicted age gradient in the stellar population inte- 
rior to LMC-4 has not been observed in the most 
recent work (Olsen et al. 1997; Braun et al. 1997; 
Dolphin & Hunter 1998); however, stochastic self- 
propagating star formation (SSPSF) is often still 
believed to play some role in the formation of 
LMC-4. An age gradient would not be produced 
if star formation within LMC-4 initiated over a 
large area at the same time, as in the bow shock- 
induced star formation scenario of de Boer et al. 
(1998). 

The models above focus mostly on the size of 
LMC-4's central cavity, and less so on the sub- 
structure of the interior stellar population. Efre- 
mov & Elmegreen (1998) tackled the problem of 
the substructure by considering the formation of 
the Quadrant and Sextant stellar arcs. They 
found evidence for a concentration of 30 Myr-old 
stars close to the center of the Quadrant arc; if 
these stars represent the remnants of a stellar clus- 
ter, they could provide the necessary pressure to 



sweep up a shell of H I gas out of which Quad- 
rant could have condensed. Efremov & Elmegreen 
(1998) thus showed that smaller-scale star forma- 
tion events have contributed, in some degree, to 
the formation of LMC-4. 

For a comprehensive understanding of how 
LMC-4 formed and has evolved, with attention 
paid to both the size of the supershell and its 
substructure, it will be necessary to survey its 
entire massive star population, measure the prop- 
erties of all of its individual H II regions, and 
study the relationship of the stars and ionized gas 
with the full distribution and kinematics of the 
neutral atomic and molecular gas. While many 
of these pieces are coming into or are already in 
place, much remains to be done, particularly in 
the areas of stellar spectroscopy and the analysis 
of large existing datasets. This work takes a small 
step in the direction of taking a comprehensive 
look at LMC-4 by focussing on the individual OB 
association LH 72, its embedded H II region DEM 
228, and its associated H I gas. Because of its iso- 
lation, LH 72 is a relatively uncomplicated object 
with which to begin the analysis. In addition, as 
remarked by Olsen et al. (1997), its location in 
LMC-4's interior may yield a special clue as to 
the formation and evolution of the surrounding 
supershell. Figure 1 displays images of LH 72 and 
DEM 228, showing their relationship to SGS-14 
and LMC-4. 

2. Observations 

Optical imaging and spectroscopy of LH 72 and 
narrow-band imaging of the embedded H II region 
DEM 228 were performed during three runs with 
the CTIO 1.5-m and 0.9-m telescopes. The pri- 
mary broad- and narrow-band imaging was car- 
ried out on the nights of November 24-26, 1998 
with the 1.5-m Cassegrain Focus imager at f/7.5 
using the SITe 2048x2048 CCD #6. The f/7.5 
secondary mirror produces a pixel scale of 0'.'44 
pixel" yielding a 15'xl5' field of view at the ex- 
pense of a somewhat erratically varying PSF. We 
observed through 3 x 3-inch Johnson UBV filters 
and 2 x 2-inch narrow-band Ha and continuum fil- 
ters. The 2-inch filters vignette the edges of the 
field, leaving a ~10'xl0' useful area. Total ob- 
ject exposure times were 1680 sec in [/, 690 sec in 
B, and 600 sec in V. We also observed Landolt 
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(1992) and Graham (1982) standards at a variety 
of airmasses each night so as to calibrate the UBV 
photometry. Calibration of the Ha images was 
achieved through 1.5-m and 0.9-m observations of 
spectrophotometric standards from Hamuy et al. 
(1994) and selected H II regions in the LMC and 
NGC 6822. 

Spectra of 12 of the brightest blue stars in 
LH 72 were taken with the 1.5-m R-C spectro- 
graph and Loral 3K detector on the nights of De- 
cember 30, 1998 and January 2-5, 1999. Using 
grating #58 in second order and a CUSO4 fil- 
ter to block contamination from first order, our 
spectra cover the wavelength range 3600-5000A, 
spanning the classical spectral classification region 
of 4000— 47OOA. Exposure times were chosen to 
achieve S/N~ 75 — 100 per resolution element. 
Comparison lamp exposures were taken at reg- 
ular intervals to provide wavelength calibration. 
Observations of the spectrophotometric standards 
EG 21 and Hiltner 600 (Hamuy et al. 1994) were 
used to provide flux calibration. 

3. The HI data 

The H I data analyzed here comes from a com- 
bination of the ATCA 21-cm survey of the Magel- 
lanic Clouds (Kim et al. 1998) and Parkes single- 
dish observations of the LMC. The mosaic of inter- 
ferometer data gives a spatial resolution of 1' and 
velocity resolution of 1.6 km s~^. The single dish 
observations provide for the zero-spacing fiux, so 
that the total mass in H I may be accurately cal- 
culated. For this project, we worked on a subset of 
the H I data cube covering a 48'x48' region cen- 
tered on 5'^23'"34!7 and -69°44'22" (2000.0) and 
spanning heliocentric velocities between 190 and 
390 km s"^ 

4. Basic reduction procedures 

AH of the CTIO frames were overscan-corrected, 
bias-subtracted, and flat-fielded using IRAF rou- 
tines with calibration data taken at the beginning 
and end of each night. For the imaging program, 
we used twilight sky flats to correct for pixel-to- 
pixel sensitivity variation. Because many of the 
flat-field and standard star exposure times were 
necessarily short, we applied a shutter shading cor- 
rection to images with exposure times less than 20 
seconds. The shutter shading map was derived by 



comparing 20-second R dome flats with R dome 
exposures during which the shutter was repeat- 
edly opened for one second and then closed. We 
also produced a bad pixel mask from short- and 
long-exposure R domes, which we applied to the 
final reduced images. 

We discovered during the later analysis that all 
of the 1.5-m CCD images had been written to the 
data acquisition computer missing the central row 
of pixels. The problem manifested itsc;lf as a 1- 
pixel shift in the positions of stars along the y-axis 
between dithered exposures. To account for the 
dropped row, we manually inserted a row of pixels 
into all of the images, while flagging the row in the 
bad pixel mask. 

For the spectral program, we used quartz lamp 
dome flats to correct for pixel sensitivity. Before 
applying the flats to the data, we corrected them 
for the shape of the quartz lamp spectrum and for 
the slit illumination function by fltting the flats 
with a cubic spline surface. Obvious cosmic rays 
were manually removed from the flat-fielded data 
by linearly interpolating over the bad pixels along 
the spatial direction. Following experimentation 
with its parameters, we used the IRAF package 
doslit to extract wavelength- and flux-calibrated 
spectra. While flux calibration may seem unnec- 
essary for spectral classification work, we found it 
easier to fit the continuum after the flux calibra- 
tion was applied, as the spectral response curve is 
wiggly compared to the smooth continuum of the 
tail of the blackbody spectrum. Our final spectral 
analysis was performed on continuum flattened 
spectra, boxcar smoothed with a 3-pixel-wide ker- 
nel. For the stars with multiple observations, we 
combined the continuum-flattened extracted spec- 
tra before smoothing. 

5. Analysis of the stellar content 
5.1. Spectral classification 

The stellar spectra were classifled indepen- 
dently by two of the authors (KO and JB). Follow- 
ing the criteria of Walborn & Fitzpatrick (1990), 
we identifled flrst the spectral type and next the 
luminosity class. Briefly, for the O stars ob- 
served here, comparison of the strengths of Hell 
4541A to Hel 4471A and HeH 4200A to Hel-t-II 
4026A provide the temperature; the strengths of 
SilV 4089A and Hell 4686A are the main lumi- 
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nosity indicators. For the early B stars, the tem- 
perature decision depends mainly on the relative 
strengths of Sill 4552A and SilV 4089A, while 
luminosity is decided through SilV 4089A, Silll 
4552A, or the strength of Mgll 4480A. In all cases, 
our identifications agreed to within two, and usu- 
ally one, spectral subtypes and one luminosity 
class. 

Figures 2a&b show the continuum-subtracted 
spectra of the 12 stars along with identified lines. 
Three of the stars, LIO, L39, and S132, show nar- 
row emission lines indicative of circumstellar disks. 
Table 1 lists our identifications. Where appropri- 
ate, the spectral type represents the average of our 
independent judgings. 

5.2. Stellar photometry 

UBV photometry of LH 72 was performed 
with the DAOPHOT and ALLSTAR suite of pro- 
grams (Stetson 1987). Aperture and point spread 
function (PSF)-fitting photometry was performed 
on individual frames, with results from multiple 
franic;s averagcxl at the end. We first performed 
aperture photometry, using a 7" radius for the 
largest aperture and 1'.'75 for the smallest. Us- 
ing DAOGROW, wc produced a growth curve 
to adjust the magnitudes measured within the 
smaller apertures so as to agree with the largest. 
DAOGROW produced a single list of magnitudes 
for each star, using the largest possible aperture 
with photometric error a < 0.02 magnitudes, 
with the magnitude extrapolated to an aperture 
of 14" radius. 

We next performed PSF photometry with ALL- 
STAR. Approximately 40 PSF stars were chosen 
automatically by DAOPHOT, the list of which 
was then edited by hand to select only uncrowded 
stars. A PSF was derived for each image by fit- 
ting a Moffat function to the stars and employing 
a quadratically varying empirical lookup table to 
correct for deviations from the analytic function. 
Aperture corrections were derived from the pre- 
viously measured aperture magnitudes, extrapo- 
lated to 14" radius. Finally, we merged all of the 
aperture and PSF magnitudes through Stetson's 
COLLECT, which produced a weighted average 
magnitude for each star in each filter, using the 
PSF magnitudes whenever those photometric er- 
rors were found to be smaller. 



The instrumental magnitudes were placed on 
the standard UBV system through photometry 
of Landolt (1992) standard fields SA 92, SA 98, 
SA 101, and T Phc and Graham's (1982) E2 and 
E3 fields. The standard star photometry was per- 
formed as for the object frames, including PSF 
photometry, as the faintest stars were underex- 
posed. We derived transformations through least- 
squares fitting to the equations: 

u = U + Ao + Ai{U-B) + A2{B-Vf + A3X + A4t 

b = B + Bo + Bi{B~V) + B3X + Bit 

v = V + Co + Cr{B-V) + CzX + dt 

where lower-case letters represent instrumental 
magnitudes, upper-case letters standard magni- 
tudes, X is the airmass, and t the time of ob- 
servation in hours from the middle of the night. 
The second-order color term A2 for U is necessary 
because of the steep response curve of the CCD in 
the ultraviolet. While usually small, variation of 
the zero point with time may make a difference at 
the level of a few hundredths of a magnitude. 

After deriving the color coefficients Ai,A2,Bi, 
and Ci independently for each night, we took their 
average and used them to rederive the remaining 
nightly coefficients. We found the average scatter 
around the transformation relations to be 0™035 
in U, 0"015 in B, and 0™008 in V. Table 2 lists 
the coefficients for the three photometric nights. 

To place the LH 72 photometry on the UBV 
system, we first transformed the photometry de- 
rived from frames taken on 11/25/98 using the re- 
lation computed for that night. We then selected 
25 isolated stars from the LH 72 frame and used 
their computed UBV magnitudes to transform 
the photometry taken on the one non-photometric 
night. Finally, all of the standardized magnitudes 
were averaged to produce a master photometric 
table for LH 72. 

Figure 3 shows the comparison of our pho- 
tometry with that from the Magellanic Clouds 
Photometric Survey (MCPS; Zaritsky, Harris, 
& Thompson 1997), kindly provided by Den- 
nis Zaritsky. Our photometry agrees within 3a 
of that of the MCPS: U - J/mcps = 0.005 ± 
0.010, B-Bmcps = 0.017±0.009,andF-VMCPS = 
0.033 ± 0.009. We also find good agreement 
with the photographic photometry of Lucke 
(1972): V - Vl72 = 0.022 ± 0.007 (39 stars), 
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{B-V)-{B-V)l72 = 0.01 ±0.02 (39 stars), and 
{U -B)-{U - B)l72 = -0.05 ± 0.03 (5 stars). 

An additional check on the photometry is to 
compare the intrinsic {U — B)o color expected 
from the spectral classification to that derived 
from the photometry after correcting for redden- 
ing. Fignrc 4 shows the comparison; we find 
([/ - S)o,spoc -{U- 5)o,phot = -0.03 ±0.01, sug- 
gesting good agreement between the photometry 
and spectroscopy. There is, however, a possible 
discrepancy at the bluest and reddest ends of the 
scale. 

5.3. Analysis of stellar photometry 

5.3.1. LH 72 membership 

Figure 5a shows the raw B — V,V color- 
magnitude diagram (CMD) produced from our 
photometry. As is typical of ground-based com- 
posite LMC CMDs, we see a broad main sequence 
down to y ~ 21, an extended giant branch, and a 
prominent red clump. The presence of an OB asso- 
ciation is clearly evidenced by the well-populated 
upper main sequence; however, the membership 
of LH 72 is confused by the background field pop- 
ulation. To restrict our analysis to probable as- 
sociation members, we selected stars lying within 
the boundary defined by Lucke (1972). Defining 
the remaining stars as representatives of the back- 
ground population, we subtracted the background 
CMD from the CMD of stars lying within LH 72 's 
boundary, as done e.g. by Olsen (1999). Figure 5b 
shows the result of the subtraction procedure. We 
sec that most of the lower main sequence, giant 
branch, and red clump have disappeared, while 
most of the upper main sequence remains. At 
^ < 16, ^80% of the stars are probable members 
of LH 72; between 16 < V < 17, the contam- 
ination rate is ~50%. For the purposes of our 
analysis, therefore, we choose as LH 72 members 
only those stars with F < 16 in the background- 
subtracted CMD. 

5.3.2. Identification of Be stars 

As discussed by e.g. Maeder, Grebel, & Mer- 
milliod (1999), the presence of emission lines in 
the spectra of late O and early B stars (the Be star 
phenomenon) is indicative of circumstellar disks. 
Fast rotation and the presence of circumstellar 
material causes Be stars to generally be redder 



and cooler than their non-Be counterparts, a fact 
which unaccounted for may be misinterpreted as 
an age difference. 

As done by e.g. Grebel (1997), we identi- 
fied Be star candidates photometrically by look- 
ing for the presence of excess Ha emission. Ha 
magnitudes were measured from the calibrated, 
continuum-subtracted Ha image of DEM 228 (dis- 
cussed in section 6) through aperture photometry 
with DAOPHOT. Figure 6 plots the Ha magni- 
tudes versus B — V. The spectroscopically iden- 
tified emission-line stars LIO and S132 show an 
Ha excess that separates them from the bulk of 
the spectroscopically studied stars. Thus, we la- 
bel those stars with Ha magnitudes < 1.6 and 
B — V < 0.5 as possible Be stars. Figure 7 shows 
that the stars selected according to this definition 
are indeed redder than the bulk of the main se- 
quence. However, we do expect some contamina- 
tion from non-Be stars, as SI 28, which does not 
contain obvious Balmer emission, is included as a 
Be candidate by our definition. 

Because we apply a magnitude cutoff of y < 16 
for certain LH 72 membership, corresponding to 
My ^ —2.5 and spectral type ~B1 on the main 
sequence, we cannot compute the Be star frac- 
tion of LH 72 over the full range of possible spec- 
tral types. We note, however, that the fraction is 
at least 10% from the number of Be stars found 
within the boundaries of LH 72. 

5.4. Reddening and the HR diagram 

We computed the reddenings, temperatures, 
and luminosities of individual stars through the 
transformations employed by Massey, Water- 
house, & DeGioia-Eastwood (2000); these are 
based on the calibrations of Vacca, Garmany, & 
ShuU (1996) for the O stars and a combination of 
sources, including Humphreys & McElroy (1984), 
as well as more modern work, for the B stars. For 
the stars with spectra, we applied the Vacca, Gar- 
many, & ShuU (1996) and Humphreys & McElroy 
(1984) calibrations directly; we then used the spec- 
tral type-intrinsic color relationship of Fitzgerald 
(1970) and our measured B — V colors to derive 
E{B - V). 

For those stars observed only in UBV, we 

used Massey et al.'s (2000) fits of T^ff to the 
reddening- free index Q. We picked a luminos- 
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ity class for each star by plotting Sclimidt-Kalcr's 
(1982) magnitude-color relationship on our CMD, 
assuming (m— M)lh72 = 18.4. While we found su- 
pcrgiants easy to distinguish from main sequence 
stars, separating giants from main sequence stars 
was made difficult by the differential reddening. 
Over the range of observed Q, confusing a giant 
with a main sequence star will introduce an error 
of up to ~0.05 dex in log Tes and ~0.3 magnitudes 
in Mboi- With Tea and the bolomctric correction 
known, we again used the spectral type-intrinsic 
color relationship of Fitzgerald (1970) to derive 
E{B-V). 

Figure 8 shows the E{B — V) distribution for 
those LH 72 stars with spectra or Q < —0.4. The 
distribution has a mean of E{B — V) = 0.09, me- 
dian of 0.075, standard deviation of 0.055, and a 
sharp lower bound at E{B -V) = 0.04. The typ- 
ical error on an individual value is ±0^01, indi- 
cating that the width of the reddening distribu- 
tion is produced by differential reddening within 
the LMC and not purely photometric errors. The 
typical foreground reddening of E{B — V)= 0.075 
towards the LMC (Schlegel, Finkbeiner, & Davis 
1998) agrees with our median value but not with 
the minimum. Outside the boundaries of LH 72, 
the reddening also has a lower bound E{B — V) = 
0.04, with the majority of stars having E{B—V) = 
0.07 ±0.04. 

Figure 9 shows the combined HR diagram for 
stars with spectra and stars with photometry 
only, overlaid with stellar evolutionary tracks and 
isochrones from Schaerer et al. (1993). With the 
exception of the 06V star S133, all of LH 72's stars 
are less massive than 40 Mq and appear older than 
4—5 Myr. However, as has been seen in a rnim- 
ber of LMC OB associations, LH 72 appears to 
contain a considerable spread in age: the H-R di- 
agram suggests that LH 72 contains stars as old 
as ~30 Myr. We see no evidence for the gradient 
in age that we reported in Olsen et al. (1997); we 
surmise that the gradient was an artifact of the 
small number of stars analyzed in that paper. 

There are a number of ways besides an age 
spread within the association to produce a broad- 
ened sequence in the observed H-R diagram. Star- 
to-star differences in the unknown rotation rates 
may explain the broadening. As shown by Heger 
& Langer (2000), fast rotation moves main se- 
quence stars to cooler temperatures and lower lu- 



minosities, as compared to the non-rotating main 
sequence. However, the addition of rotation to a 
stellar evolutionary track also increases the main 
sequence lifetime. The increase in lifetime is such 
as to produce a real age spread of > 15 Myr if 
the broadening in LH 72's H-R diagram were due 
entirely to stellar rotation, a number very simi- 
lar to what we claim on the basis of non-rotating 
stellar models. A second possibility is that line-of- 
sight contamination within the LMC produces the 
broadened sequence. However, because our anal- 
ysis shows that 80% of the OB stars with y < 16 
lie within the confines of LH 72's boundary, we 
consider this possibility unlikely. 

A third, more likely, possibility is that a sys- 
tematic error in the photometry is responsible for 
the broadened sequence. During our comparison 
of the intrinsic colors of the stars with spectra with 
their observed Q values, we found good agreement 
of the colors derived from spectroscopy with those 
from photometry. However, if we assume that we 
should lower the Q's of stars with only photome- 
try by as much as Of^OB, then we produce the H-R 
diagram of Figure 10. The shift in Q makes LH 
72 appear younger; the envelope of the youngest 
stars now suggests an age of ~3 Myr instead of 
5 Myr. The apparent age spread is reduced by a 
factor -^2; the oldest stars now have ages of ~20 
Myr, suggesting an age spread of ~15 Myr. Such 
an age spread is still significant compared to the 
measurement errors in Tes and Mboi- 

6. Analysis of the ionized gas 
6.1. Calibration 

The individual Ha and off-band images of DEM 
228 were first background-subtracted. We used a 
400 pixel-wide box placed in one corner of each 
image to measure the modal sky value. The obser- 
vations of spectrophotometric standards were next 
used to measure the Ha zero point and extinction 
coefficient. To correct for the different illumina- 
tion of the filter passband by the standards and 
by the H II region, we used CTIO's published fil- 
ter transmission curve and the Hamuy et al. (1994) 
published spectra to calculate the number of pho- 
tons detected at the LMC Doppler-shifted wave- 
length of Ha. The flux of LMC Ha photons de- 
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tected through observation of a standard star is: 



NHa = N 



J F{\)T{X)S{XHa)dX 

jF{X)T{X)dX 



where N is the flux of photons observed through 
the filter, F{X) is the calibrated spectrum of the 
star, T(A) is the filter transmission curve, and Ha 
has been redshifted to the mean velocity of LH 72 's 
associated B.I,v = 307 km s'^ (Kim et al. 1998). 
We used quadratic interpolation to bring Hamuy 
et al.'s spectra, which are tabulated at 50 A in- 
tervals, to the wavelength grid of T{X). We then 
measured the zero point and extinction coefficient 
by fitting for Zna and kna in the equation: 

2.51ogio^ffa = ZHa+kHaX+2.5logio j F {X)S{XH<,)dX 

where X is the airmass. We measured kna = 
0.08 ± 0.01, in good agreement with the usual R 
extinction coefficient of 0.085 at Cerro Tololo. 

After applying the calibration to the individ- 
ual Ha images, we removed cosmic rays and 
subtracted the off-band images. Cosmic rays 
were identified in both Ha and off-band images 
with the IDL Astronomy User's Library routine 
CR_REJECT, and removed by linear interpola- 
tion through the surrounding good pixels. We 
next subtracted the nearest available off-band im- 
age, shifting the image to adjust for fractional 
pixel offsets and scaling it to match the fluxes 
of point sources through the Ha filter. Finally, 
we shifted and median combined the calibrated, 
continuum-subtracted Ha images into a single im- 
age. 

6.2. Correction of Ha flux for reddening 

As shown by the UBV photometry, variable 
reddening across DEM 228 is significant. One 
way to correct for the spatially varying redden- 
ing would be to observe the Hf] emission line and 
to construct a map of the Ha/H/3 ratio. However, 
in the case where the dust is embedded in the ion- 
ized gas, scattered nebular light complicates the 
interpretation of the Ha/H/3 ratio (Mathis 1970; 
Skelton 1999). Instead, we chose to use the red- 
dening values derived for the stars themselves to 
construct a reddening map. 

We produced the reddening map through lin- 
ear interpolation of the discrete reddening val- 



ues. We first constructed a mesh of triangles con- 
necting nearby stars using Delaunay triangulation 
(through the IDL procedure TRIANGULATE), 
tiling the surface boimded by the locations of the 
observed stars. We then linearly interpolated the 
reddening map over the grid using the values at 
the vertices. Figure 11 shows our Ha image over- 
laid with smoothed contours of the reddening map. 
The map shows that the brightest regions of DEM 
228 tend to have higher reddening, but that the 
northern portion of the nebula is more heavily ob- 
scured than the southern end. The peak redden- 
ing value is E{B — V) = 0.23 while we assume 
E{B — V) = 0.04 for the background reddening. 
We used the reddening curve published by Mathis 
(1990) to correct the Ha surface flux for absorp- 
tion by dust. 

6.3. Comparison of Ha flux with stellar 
ionizing radiation 

Is the Ha emission from DEM 228 fully ac- 
counted for by LH 72 's O and B stars and those of 
the field? We estimated the luminosity of Lyman 
continuum photons Qo produced by the observed 
hot stars using the calibration of Schaerer & de 
Koter (1997). This calibration is based on theoret- 
ical spectral energy distributions of massive stars, 
calculated using models incorporating non-LTE ef- 
fects, stellar winds, and line blanketing. The pa- 
rameters affecting Qo are temperature, gravity, 
and metallicity. To derive the gravities, we es- 
timated the masses by fitting to the Z = 0.008 
Geneva tracks (Schaerer et al. 1993) and calcu- 
lated the radii through L = AirR'^aT^ff, L and T^s 
being known through our spectra and photome- 
try. Next, we identified all entries in Schaerer 
& de Koter's (1997) Table 3 with temperatures 
within the errors of the temperatures of the ob- 
served stars; from that group of entries, we chose 
the entry with log g closest to the calculated log g. 
Finally, we interpolated Qo for the metallicity of 
the LMC from the Z = 0.02 and Z = 0.004 entries. 

Assuming Case B recombination (Osterbrock 
1989) and an LH 72 distance modulus of (m — 
M)o = 18.4, we next derived the Ha fluxes that 
would be produced by the computed Lyman con- 
tinuum luminosities. To get an idea of the contri- 
bution of each star to the morphology of the H II 
region, we then identified those circular regions in 
the de-reddened Ha image which can account for 
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all of the stellar ionizing radiation, taking into ac- 
count all stars simultaneously. Figure 12 shows 
the de-reddened Ha image with the circular re- 
gions overlaid. The two most luminous stars, S128 
and S133, are marked with X's. Their Lyman 
continuum photons account for more Ha emission 
than is measured in the imago; DEM 228 thus ap- 
pears to be a density-bounded H II region. The 
remaining stars account for the bulk of the high- 
surface brightness emission. Most of this emission 
appears to derive from multiple sources; for exam- 
ple, at least 15 sources contribute to the ionization 
of the southernmost complex of H II regions. On 
the other hand, two of the bright knots appear to 
be single-star H II regions. 

Although the total reddening-corrected Ha 
emission accounts for only 60% of the Lyman con- 
tinuum photons produced, not all of the structure 
in the ionized gas is readily explained by the distri- 
bution of hot stars. In particular, the low surface 
brightness halo of Ha emission contains structure 
suggesting past star formation activity. If true, 
this implies that the area hosting star formation 
within LH 72 has shrunk somewhat and drifted 
northwards. 

7. HI structure and kinematics 

Figure 13 shows the individual channel maps of 
the HI cube covering the velocity range 272.5 < 
V < 320.3 km s'^ The eighth-note shape of DEM 
228 is mimicked by the nearby HI gas, most promi- 
nently between the velocities 300 ^ v < 310 km 
s~^. In Kim et al. (1999), we identified this gas 
as forming the western edge of a supergiant shell, 
SGS-14, an observation also made by Efremov & 
Elmegreen (1998). However, because of the appar- 
ent lack of evidence for an approaching or receding 
cap to the shell, we did not then report an expan- 
sion velocity for it. 

We now propose that the surface of SGS-14 is 
not smooth but patchy, and that the HI data in- 
deed suggest that the shell is expanding. We base 
this claim mainly on the evolution of the HI gas 
associated with DEM 228 through the sequence of 
channel maps shown in Figure 13. In the v = 300.5 
km s~^ map, an outline similar in shape to DEM 
228 is clearly seen. Following the northeastern tip 
of the outline, we notice that the tip moves pro- 
gressively towards the interior of the shell with 



increasing velocity. The first position- velocity di- 
agram shown in Figure 14 also demonstrates this 
behavior of the gas, indicating that the filament 
associated with DEM 228 is expanding with a ve- 
locity of ~ 15 km s~^. 

The channel maps showing the eastern rim of 
SGS-14 also suggest expansion. In the v = 302 km 
s~^ map, the faint trace of the "C" -shaped eastern 
shell edge can be seen. The arms of the "C" pinch 
together with increasing velocity, as would be ex- 
pected if the gas formed part of an expanding shell. 
The "C" fades into the noise at u ~ 315 - 320 km 
s^^, roughly the same velocity at which the tip of 
the filament forming the western edge of SGS-14 
disappears. The second position- velocity diagram 
shown in Figure 14 implies an expansion velocity 
of ~10 km s""'^. 

The u ~ 274 — 284 km s~^ maps contain what 
may be a piece of the approaching portion of 
SGS-14. The elongated cloud of gas lies roughly 
halfway between the geometric center of SGS-14 
and its incomplete southern boundary. The mean 
velocity of the cloud, v ^ 280 km s"-*^, suggests an 
expansion velocity of ~20 km s~^, as shown in the 
third panel of Figure 14. 

The mass of HI that is involved in the 10-20 km 
s~^ expansion can be calculated from the observed 
column density of the gas. The lower left panel 
of Figure 1 shows the assumed physical bound- 
ary of SGS-14. Integrating over this area and 
the velocity range (280,320) km s~^, we find a 
mass 

MsGS-14 = 1-8 X 10^ Mq, with 3.3xlO''M0 
purely in the gas associated with LH 72. The ki- 
netic energy carried by the H I of SGS-14 is thus 
^ 4x 10"'''^ ergs, or that produced by ~4 supernovae 
if one assumes that 10% of the explosion energy 
(^ 10^^ ergs) goes into moving the surrounding 
ISM (Tenorio-Tagle 1991). 

Assuming a distance modulus (m — M) = 18.4, 
the radius of SGS-14 is ~240 pc. Thus, if the 
gas in SGS-14 were uniformly distributed over the 
sphere, it would have an average density po ~ 0.1 
cm~^. If SGS-14 expanded with constant velocity 
of ~15 km s~^ over its entire lifetime, its age is 
~16 Myr. 

8. Discussion and conclusion 

Having examined the stellar content of LH 72, 
studied the properties of the ionized gas of the 
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embedded H II region DEM 228, and derived the 
mass and kinematics of the surrounding H I super- 
shell SGS-14, a plausible story emerges: SGS-14 is 
an expanding supcrshcU which has cleared the NE 
quadrant of LMC-4 of gas and triggered the for- 
mation of the OB association LH 72 behind its 
shock front. 

What is the evidence? First, the ages of the 
oldest stars in LH 72, ~ 15 — 30 Myr, agree closely 
with the expansion age of SGS-14, ^15 Myr. Since 
the first triggering by the passing supershell, star 
formation could then have continued in LH 72 un- 
til < 5 Myr ago, although without the clear N-S 
spatial pattern reported by Olsen et al. (1997). 
Second, the interior of SGS-14 appears empty of 
gas; the H I maps show that the 21-cm emission 
is at the level of the noise. Moreover, DEM 228 is 
density-bounded, allowing ionizing photons from 
the UV-bright stars S128 and S133 to escape into 
the interior of SGS-14. Third, overlaying our Ha 
image with contours of H I column density show 
that while DEM 228 and LH 72 he close to the rim 
of SGS-14, they are displaced towards its interior. 
Figure 15 shows the relationship. The plate solu- 
tion for the Ha image was calculated using stars in 
LH 72 with coordinates from SIMBAD. The calcu- 
lated positions of four of the brightest stars agree 
to within < 1" with those listed in the Tycho-2 
catalog (H0g et al. 2000); the astrometric solution 
should thus be excellent. The HI astrometric so- 
lution is accurate to << 1'. The figure shows that 
the peak Ha emission is offset by 1-2 arminutes 
(15-30 pc) from the peak HI column density, to- 
wards the interior of SGS-14. 

This picture of the formation of SGS-14 and LH 
72 has an interesting consequence for our under- 
standing of the supergiant shell LMC-4. LMC-4 
has often been thought to have formed as a unit 
(e.g. Westerlund & Matliewson 1966; Dopita et al. 
1985; Efremov et al. 1998). However, this work 
shows that one quadrant of LMC-4 consists of an 
independently expanding shell, which may have 
triggered the formation of at least one of LMC- 
4's OB associations, LH 72, along its rim. Thus, 
a possible formation scenario of LMC-4 is one in 
which many smaller shells, some triggering star 
formation from the swept-up gas, have overlapped 
to produce the ring of H I and H II regions we see 
today. In this scenario, the energy required to pro- 
duce LMC-4 is greatly reduced from that needed 



if it formed as a single unit-the reason being that 
once a supershell blows out, it becomes difficult to 
impart additional energy to expansion within the 
disk plane (e.g. Mac Low & McCray 1988). LMC- 
4, if it formed as a unit, must have blown out long 
ago, as its diameter of ~1.4 kpc is many times the 
scale height of the LMC's H I disk of 180 pc (Kim 
et al. 1999). By contrast, with a radius of ~ 240 
pc, SGS-14 should be only in the early stages of 
blowout. 

The scenario proposed here naturally merges 
with that of Efremov & Elmegreen (1998), in 
which Quadrant and Sextant formed from swept- 
up shells which have since merged with SGS-11. 
Indeed, LH 72 may be a younger, less massive 
version of Quadrant and Sextant. Like those two 
structures, it is arc-shaped, with an approximate 
center of curvature near the center of SGS-14. In 
^15 Myr, by the time that the last of the cur- 
rent generation of ionizing stars have disappeared, 
SGS-14 will have approximately doubled in size, 
assuming constant expansion velocity. Assuming 
that LH 72 also grows only in the azimuthal di- 
rection, it will still have an overdensity of stars at 
the tip of the main sequence compared to the field, 
by a factor ~4. LH 72 may then be classified as 
a stellar arc, shorter and less tightly curved than 
the current Quadrant. 

We are grateful to the NOAO Time Allocation 
Committee for the time awarded to this project, 
and to the CTIO mountain staff for their excellent 
support. We thank Bruce Elmegreen for enlight- 
ening discussions, and Paul Hodge for the ideas 
that led to this project. We thank the anonymous 
referee for prompt and helpful comments. 
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Fig. 1.— 

A montage showing the relationships of LH 72 
and DEM 228 to SGS-14, LMC-4, and the LMC. 
Proceeding counter-clockwise from top left, the 
images are: Boyden Observatory photograph of 
the LMC, courtesy of Paul Hodge and Harlow 
Shapley; same photograph with closc-up of region 
containing LMC-4; section of H I map (Kim et al. 
1998) showing the location of SGS-14; true-color 
UBV combined image of OB association LH 72 
taken with the CTIO 1.5-m telescope; CTIO 1.5- 
m Ha image of the H II region DEM 228. 

Fig. 2. — (a) continuum- flattened, wavelength- 
calibrated spectra of 9 O and B stars in LH 72. 
All spectra are plotted on the same arbitrary flux 
scale. The locations of a number of typical absorp- 
tion lines are indicated by vertical lines. (6) as in 
a, but showing stars containing Balmer emission 
lines. 

Fig. 3.- 

Comparison of our UBV photometry with that 
from the Magellanic Clouds Photometric Survey 
(Zaritsky, Harris, & Thompson 1997). Points with 
positive values along the y axes indicate stars for 
which our photometry is fainter than that of the 
MCPS. The solid lines are medians through the 
points, while the dotted lines are la deviations of 
the means. 

Fig. 4.— 

Comparison oi U — B colors derived from our 
photometry, after correction for reddening as de- 
scribed in the text, with the intrinsic U — B colors 
derived from our spectral classifications. For most 
of the stars, the colors derived from the spectra 
are bluer than those found from the photometry, 
although the median difference is only ~0.03 mag- 
nitudes. 

Fig. 5. 

(a) B — V,V color-magnitude diagram of LH 
72 and its surrounding field. Error bars, shown 
function of V magnitude, are the median la 
photometric errors reported by DAOPHOT. (6) 
Color-magnitude diagram containing stars found 
only within the borders of LH 72, as defined by 
Lucke (1972). Small points are stars removed by 
the subtraction procedure, while large circles are 
those that remain (see text for details). 



Fig. 6. - 

Ha magnitudes, defined as — 2.5 log(///a) -|- 25, 
plotted versus B — V; fua is measured in units 
of 10""'^^ ergs cm~^ s~^. Stars with spectra are 
labelled with filled circles. Those stars with Ha 
fluxes greater than or equal to those of LH72-1, 
LIO, and S128 and with B—V < 0.5 are considered 
to be candidate Be stars. 

Fig. 7.— 

LH 72 color-magnitude diagram showing the lo- 
cations of the candidate Be stars (fllled circles) se- 
lected through Figure 6. The Be candidates have 
red colors compared to the main sequence. 

Fig. 8.— 

Histogram of reddening values E[B — V) for 
stars within LH 72. The distribution has a me- 
dian of E{B -V) = 0.075, mean of 0.09, standard 
deviation of 0.055, and a sharp lower bound at 
{B — V) = 0.04. The typical error on an individ- 
ual measurement is 0™01. 

Fig. 9.— 

Hertzprung-Russell diagram of stars in LH 72, 
excluding candidate Be stars. Black circles indi- 
cate stars with both spectroscopy and photometry, 
gray circles stars with photometry only. The grey 
lines are evolutionary tracks from Schaerer et al. 
(1993), with initial masses labelled to the left of 
the tracks. The black lines are isochrones with 
ages of 1, 2, 4, 8, 16, and 32 Myr. 

Fig. 10.— 

H-R diagram produced if the observed Q val- 
ues are lowered by 0.05 magnitudes. Evolutionary 
tracks and isochrones are as in Figure 9. 

Fig. 11.— 

Ha image overlaid with smoothed contours of 
the reddening map, which was derived by interpo- 
lating over stellar E{B — V) values as described in 
the text. The contours are labelled by E{B — V); 
lines perpendicular to the contours indicate the 
direction of the local downhill gradient. 

Fig. 12.— 

Ha image, stretched by the square root of the 
intensity, overlaid with circles indicating the re- 
gions which account for all of the Lyman contin- 
uum emission from the enclosed hot stars (see text 
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for an explanation of the calculation). The stars 
S128 and S133 are marked by X's-their Lyman 
continuum contributions exceed the available Ha 
flux. The observed Ha emission accounts for only 
60% of the available Lyman continuum photons, 
indicating that DEM 228 is density-bounded. 

Fig. 13.— 

Individual H I channel maps of the region sur- 

roimding SGS-14, spaced by units of the velocity 
resolution. Velocities are in km s~^; the spatial 
scale is indicated in the first channel. The gas as- 
sociated with LH 72 /DEM 228 is most prominent 
in the velocity range 295.5 < v < 313.7 km s~^. 

Fig. 14.— 

Position-velocity maps at three positions. In 
each panel, the H I map of SGS-14 is shown, 
summed over a range of velocities. The white lines 
simulate a slit the velocity structure at that slit 
position is shown either above or to the right of 
the spatial map. Velocity labels are in km s~^. 

Fig. 15.— 

Ha image, stretched by the logarithm of the in- 
tensity, overlaid with contours of H I column den- 
sity in units of 10^^ atoms cm~^. North is up and 
east is to the left. 



13 



Table 1 

Spectral and Luminosity Classes of Stars in LH 72 



Star 


R.A. (2000.0) 


Dec. 


Sp. type 


Lum. class 


Comment 


S133 


5'' 32"' 29^:2 


-66° 


28' 


00'.' 2 


06 


V 




LH72-1 


5'' 32"' 27"8 


-66° 


28' 


20'.' 4 


08.5 


V 




L29 




-66° 


28' 


02" 9 


09 


III 




L41 


5^ 32" 18!4 


-66° 


25' 


51^1 


09.5 


V 




S128 


5^ 31m 53?4 


-66° 


31' 


14;'3 


09.7 


lb 




L22 


5^ 32" 18! 1 


-66° 


28' 


47" 2 


BO 


V 




L49 


5^ 32" 13!0 


-66° 


25' 


03;'0 


BO. 2 


V 




L51 




-66° 


23' 


58^' 2 


B0.5 


III 




S135 


5^ 32" 33!2 


-66° 


26' 


12" 3 


B3 


la-b 




LIO 


5^ 32" 33?7 


-66° 


27' 


lO^'O 


08 


V 


emission lines 


L39 


gh 32m Qgsg 


-66° 


26' 


08" 4 


B0.5 


V 


emission lines 


S132 


S*" 32" 18!8 


-66° 


24' 


12;'0 


B8 


la 


emission lines 



Table 2 

Photometric transformation coefficients 



Date 


^0 


Ai 


A2 


A3 


Ai 


11/25/98 

11/26 

11/27 


4.249(7) 
4.323(5) 
4.278(5) 


-0.167 
-0.167 
-0.167 


0.136 
0.136 

0.136 


0.374(60) 
0.555(42) 

0.485(25) 


0.0 
-0.009(2) 
0.0 




Bo 


Bi 


B3 


B4 




11/25/98 

11/26 

11/27 


2.414(3) 

2.426(2) 
2.391(2) 


0.112 

0.112 
0.112 


0.287(30) 

0.302(16) 
0.260(10) 


0.0 

-0.003(1) 
0.0 






Co 


Ci 


C3 


C4 




11/25/98 

11/26 

11/27 


2.084(2) 
2.102(2) 
2.091(1) 


-0.029 
-0.029 
-0.029 


0.149(18) 
0.139(14) 
0.140(5) 


0.0 
-0.005(1) 
0.0 
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